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Abstract 

We are given a set of weighted unit disks and a set of points in Euclidean plane. The 
minimum weight unit disk cover (WUDC) problem asks for a subset of disks of minimum total 
weight that covers all given points. WUDC is one of the geometric set cover problems, which 
have been studied extensively for the past two decades (for many different geometric range 
spaces, such as (unit) disks, halfspaces, rectangles, triangles). It is known that the unweighted 
WUDC problem is NP-hard and admits a polynomial-time approximation scheme (PTAS). For 
the weighted WUDC problem, several constant approximations have been developed. However, 
whether the problem admits a PTAS has been an open question. In this paper, we answer this 
question affirmatively by presenting the first PTAS for WUDC. Our result implies the first PTAS 
for the minimum weight dominating set problem in unit disk graphs. Combining with existing 
ideas, our result can also be used to obtain the first PTAS for the maximum lifetime coverage 
problem and an improved 4.475-approximation for the connected dominating set problem in 
unit disk graphs. 
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1 Introduction 


The set cover problem is a central problem in theoretical computer science and combinatorial 
optimization. A set cover consists of a ground set U and collection S of subsets of U. Each set 5 G 5 
has a non-negative weight ws- The goal is to find a subcollection C C 5 of minimum total weight 
such that y C covers all elements of U. The approximibility of the general SC problem is rather well 
understood: it is well known that the greedy algorithm is an //^-approximation (i/„ = 1 /^) 

and obtaining a (1 — e) In n-approximation for any constant e > 0 is NP-hard [15,23]. In the 
geometric set cover problem, U is a set of points in some Euclidean space and S consists 
of geometric objects (e.g., disks, squares, triangles). In the geometric setting, we can hope for 
better-than-logarithmic approximations due to the special structure of S. Most geometric set cover 
problems are still NP-hard, even for the very simple classes of objects such as unit disks [11,29] 
(see [9,27] for more examples and exceptions). Approximation algorithms for geometric set cover 
have been studied extensively for the past two decades, not only because of the importance of the 
problem per se, but also its rich connections to other important notions and problems, such as VC- 
dimension [5,12,22], e-net, union complexity [10,39,40], planar separators [25,35], even machine 
scheduling problems [3]. 

In this work, we study the geometric set cover problem with one of the simplest class of objects 
— unit disks. The formal definition of our problem is as follows: 

Definition 1. Weighted Unit Disk Cover (WUDC): Given a set D = {Di,..., D^} of n unit disks 
and a set V = {Pi, ..., Pm} of m points in Euclidean plane M^. Each disk Dj has a weight rc(Dj). 
Our goal is to choose a subset of disks to cover all points in V, and the total weight of the chosen 
disks is minimized. 

WUDC is the general version of the following minimum weight dominating set problem in unit 
disk graphs (UDG). 

Definition 2. Minimum Weight Dominating Set (MWDS) in UDG: We are given a unit disk graph 
G{V, E), where V is a set of weighted points in and {u, v) £ E iff ||u — n|| < 1 for any u,v G V. 
A dominating set S is a subset of V such that for any node v ^ S, there is some u £ S with 
{u,v) £ E. The goal of the minimum weight dominating set problem is to find a dominating set 
with the minimum total weight. 

Given a dominating set instance with point set V, we create a WUDC instance by placing, for 
each point in v £ V, a point (to be covered) co-located with v and a unit disk centered at v, with 
weight equal to the weight of v. Thus, WUDC is a general version of problem MWDS in UDG. In 
this paper, we only state our algorithms and results in the context of WUDC. 

1.1 Previous Results and Our Contribution 

We first recall that a polynomial time approximation scheme (PTAS) for a minimization problem 
is an algorithm A that takes an input instance, a constant e > 0, returns a solution SOL such that 
SOL < (1 -|- e)OPT, where OPT is the optimal value, and the running time of A is polynomially in 
the size of the input for any fixed constant e. 

WUDC is NP-hard, even for the unweighted version (i.e., rc(Dj) = 1) [11]. For unweighted 
dominating set in unit disk graphs. Hunt et al. [31] obtained the first PTAS in unit disk graphs. 
For the more general disk graphs, based on the connection between geometric set cover problem 
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and e-nets, developed in [5,12,22], and the existence of e-net of size 0(l/e) for halfspaces in [36] 
(see also [26]), it is possible to achieve a constant factor approximation. As estimated in [35], these 
constants are at best 20 (A recent result [6] shows that the constant is at most 13). Moreover, 
there exists a PTAS for nnweighted disk cover and minimum dominating set via the local search 
technique [25,35]. 

For the general weighted WUDC problem, the story is longer. Ambiihl et al. [2] obtained the 
first approximation for WUDC with a concrete constant 72. Applying the shifting techniqne of [28], 
Huang et al. [30] obtained a (6-|-e)-approximation algorithm for WUDC. The approximation factor 
was later improved to (5 -|- e) [13], and to (4 -|- e) by several groups [14,20,42]. The current best 
ratio is 3.63. ^ Besides, the quasi-uniform sampling method [10,40] provides another approach to 
achieve a constant factor approximation for WUDC (even in disk graphs). However, the constant 
depends on several other constants from rounding LPs and the size of the union complexity. Very 
recently, based on the separator framework of Adamaszek and Wiese [1], Mustafa and Raman [34] 
obtained a QPTAS (Quasi-polynomial time approximation scheme) for weighted disks in (in 
fact, weighted halfspaces in M^), thus ruling out the APX-hardness of WUDC. 

Another closely related work is by Erlebach and van Leeuwen [21], who obtained a PTAS for 
set cover on weighted unit sqnares, which is the first PTAS for weighted geometric set cover on 
any planar objects (except those poly-time solvable cases [9,27]). Although it may seem that their 
result is quite close to a PTAS for weighted WUDC, as admitted in their paper, their technique is 
insufficient for handling unit disks and “completely different insight is required”. 

In light of all the aforementioned results, it seems that we should expect a PTAS for WUDC, 
but it remains to be an open question (explicitly mentioned as an open problem in a number of 
previous papers, e.g., [2,17,19-21,38] ). Our main contribution in this paper is to settle this 
question affirmatively by presenting the first PTAS for WUDC. 

Theorem 3. There is a polynomial time approximation scheme for the WUDC problem. The 
running time is ^. 

Because WUDC is more general than MWDS, we immediately have the following corollary. 

Corollary 4. There is a polynomial time approximation scheme for the minimum weight dominat¬ 
ing set problem in unit disk graphs. 

We note that the running time jg nearly optimal in light of the negative result by 

Marx [33], who showed that an EPTAS (i.e., Efficient PTAS, with running time /(l/e)poly(n) ) 
even for the unweighted dominating set in UDG would contradict the exponential time hypothesis. 

Einally, in Section 9, we show that our PTAS for WUDC can be used to obtain improved 
approximation algorithms for two important problems in wireless sensor networks, the connected 
dominating set problem and the maximum lifetime coverage problem in UDG. 

2 Our Approach - A High Level Overview 

By the standard shifting technique [16], it suffices to provide a PTAS for WUDC when all disks 
lies in a square of constant size (we call it a block, and the constant depends on 1/e). This idea is 
formalized in Huang et al. [30], as follows. 

^ The algorithm can be found in Du and Wan [17], who attributed the result to a manuscript by Willson et al. 
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Lemma 5 (Huang ta al. [30]). Suppose there exists a p-approximation for WUDC in a fixed L x L 
block, with running time f{n,L). Then there exists a {p + 0(1/L))-approximation with running 
time 0{L-n- f{n, L)) for WUDC. In particular, setting L = 1/e, there exists a {p+e)-approximation 
for WUDC, with running time O (^ • n • f{n, ^)). 

In fact, almost all previous constant factor approximation algorithms for WUDC were obtained 
by developing constant approximations for a single block of a constant size (which is the main 
difficulty ^). The main contribution of the paper is to improve on the previous work [2,13,20,30] 
for a single block, as in the following lemma. 

Lemma 6. There exists a PTAS for WUDC in a fixed block of size Lx L for L = 1/e. The running 
time of the PTAS is 

From now on, the approximation error guarantee e > 0 is a fixed constant. Whenever we say 
a quantity is a constant, the constant may depend on e. We use OPT to represent the optimal 
solution (and the optimal value) in this block. We use capital letters A, B,C,.. . to denote points, 
and small letters a,b,c,.. . to denote arcs. For two points A and B, we use \ AB\ to denote the line 
segment connecting A and B (and its length). We use Dj to denote a disk and Di to denote its 
center. For a point A and a real r > 0, let D(H, r) be the disk centered at A with radius r. For a 
disk Dj, we use dDi to denote its boundary. We call a segment of 5Dj an arc. 

First, we guess that whether OPT contains more than C disks or not for some constant C. If 
OPT contains no more than C disks, we enumerate all possible combinations and choose the one 
which covers all points and has the minimum weight. This takes O ~ 0{n^) time, 

which is polynomial. 

The more challenging case is whether OPT contains more than C disks. In this case, we guess 
(i.e., enumerate all possibilities) the set Q of the C most expensive disks in OPT. There are at 
most a polynomial number (i.e., 0{n^)) possible guesses. Suppose our guess is correct. Then, we 
delete all disks in Q and all points that are covered by Q. Let Dj (with weight wt) be the cheapest 
disk in Q. We can see that OPT > Cwt. Moreover, we can also safely ignore all disks with weight 
larger than wt (assuming that our guess is correct). Now, our task is to cover the remaining points 
with the remaining disks, each having weight at most wt. We use V = D\Q and V = V \ V{Q) 
to denote the set of the remaining disks and the set of remaining points respectively, where 'P(^) 
denote the set of points covered by some disk in Q. 

Next, we carefully choose to include in our solution a set % T D' oi at most eC disks. The 
purpose of TL is to break the whole instance into many (still a constant) small pieces (substructures), 
such that each substructure can be solved optimally, via dynamic programming. ^ One difficulty 
is that the substructures are not independent and may interact with each other (i.e., a disk may 
appear in more than one substructure). In order to apply the dynamic programming technique to 
all substructures simultaneously, we have to ensure the orders of the disks in different substructures 
are consistent with each other. Choosing % to ensure a globally consistent order of disks is in fact 
the main technical challenge of the paper. 

Suppose we have a set TL which suits our need (i.e., the remaining instance iV' \ \ V{Tl)) 

can be solved optimally in polynomial time by dynamic programming). Let S be the optimal 

^ For the unweighted dominating set problem in a single block, it is easy to see the optimal number of disks is 
bounded by a constant, which implies that we can compute the optimum in poly-time. However, for the weighted 
dominating set problem or WUDC the optimal solution in a single block may consist of 0{n) disks. 

^ An individual substructure can be solved using a dynamic program similar to [2,27]. 
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solution of the remaining instance. Our final solution is SOL = Q S. First, we can see that 

w{S) < n;(OPT -Q-U)^ OPT - w{g), 

since OPT — Q — %\s a, feasible solution for the instance iV' \ T-L,V' \ Hence, we have that 

SOL = w{g) + w{n) + n;(5) < OPT + eCwt < (1 + e)OPT, 

where the 2nd to last inequality holds because I'Hl < eC, and the last inequality uses the fact that 
OPT > w{g) > Cwt. 

Constructing %•. Now, we provide a high level sketch for how to construct % C T)'. First, we 
partition the block into small squares with side length /r = 0(e) such that any disk centered in 
a square can cover the whole square and the disks in the same square are close enough. Let the 
set of small squares be H = {rij}i<i,j<K where K = Lj^. For a small square F, let Dsj, G F and 
Dtp G F be the furthest pair of disks (i.e., iD^r-^trl maximized). We include the pair D^p and 
Dtp in Ti, for every small square F G H, and call the pair the square gadget for F. See Figure 1 for 
an example. We only need to focus on covering the remaining points in the uncovered region \]('H). 

We consider all disks in a small square F. The uncovered portion of those disks defines two 
disjoint connected regions (See right hand side of Figure 1, the two shaded regions). We call such 
a region, together with all relevant arcs, a substructure (formal definition in Section 4). In fact, we 
can solve the disk covering problem for a single substructure optimally using dynamic programming 
(which is similar to the dynamic program in [2,27]). It appears that we are almost done, since 
(“intuitively”) all square gadgets have already covered much area of the entire block, and we 
should be able to use similar dynamic program to handle all such substructures as well. However, 
the situation is more complicated (than we initially expected) since the arcs are dependent. See 
Figure 2 for a “not-so-complicated” example. Firstly, there may exist two arcs (called sibling 
arcs) which belong to the same disk when the disk is centered in the core-central area, as show in 
Figure 1). The dynamic program has to make decisions for two sibling arcs, which belong to two 
different substructures (called R(emotely)-correlated substructures), together. Second, in order to 
carry out dynamic program, we need a suitable order of all arcs. To ensure such an order exists, 
we need all substructures interact with each other “nicely”. 

In particular, besides all square gadgets, we need to add into H a constant number of extra 
disks. This is done by a series of “cut” operations. A cut can either break a cycle, or break one 
substructure into two substructures. To capture how substructures interact, we define an auxiliary 
graph, call substructure relation graph S, in which each substructure is a node. The aforementioned 
R-correlations define a set of blue edges, and geometrically overlapping relations define a set of red 
edges. Though the cut operations, we can make blue edges form a matching, and red edges also 
form a matching, and S acyclic (we call S an acyclic 2-matching). The special structure of S 
allows us to define an ordering of all arcs easily. Together with some other simple properties, we 
can generalize the dynamic program from one substructure to all substructures simultaneously. 

3 Square Gadgets 

We discuss the structure of a square gadget Gg(F) associated with the small square F. Recall that 
the square gadget Gg(F) = D^ U Dt, where D^ and Dj are the furthest pair of disks in F. We can 
see that for any disk Dj in F, there are either one or two arcs of dDi which are not covered by 
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Figure 1: A square gadget. and are the furthest pair of disks in square F whose centers are 
Dg and Dt- On the left hand side, the blue region is the central area G1 = D{Ds,rst) H D{Dt,rst), 
where Vgt = \DsDt\. The brown region is the core-central area Co = D(P, 1) n D((5,1). On the 
right hand side, the green area is the active regions, defined as (Uiec O 0^)) n H~^ and 

(Uec„D.-(D.UDt))ni/-. 


Gg(r). Without loss of generality, assume that DgDt is horizontal. The line DgDt divides the 
whole plane into two half-planes which are denoted by (the upper half-plane) and H~ (the 
lower half-plane). and 50^ intersect at two points P and Q. We need a few dehnitions which 
are useful throughout the paper. Figure 1 shows an example of a square gadget. 


1. (Central Area and Core-Central Area) Define the central area of Gg(r) as the intersection 
of the two disks D{Ds,rst) and D{Dt,rst) in the square F, where rgt = \DsDt\. We use C to 
denote it. Since and are the furthest pair, we can see that every other disk in F is 
centered in the central area C. 

We define the core-eentral area of Gg(r) is the intersection of two unit disks centered at P, Q 
respectively. Essentially, any unit disk centered in the core-central area has four intersections 
with the boundary of gadget. Let us denote the area by Co- 

2. (Active Region) Consider the regions 



-(D.UDi) riH-' 


and 



-(D.UDt) nH~ 


We call each of them an active region associated with square F. An active region can be 
covered by disks centered in the core-central area. We use Ar to denote an active region. 


4 Substructures 

Initially, Ti includes all square gadgets. In Section 7, we will include in 77 a constant number of extra 
disks. For a set S of disk, we use M(5) to denote the region covered by disks in S (i.e., Uoie^Dj). 
Assuming a fixed Ti, we now describe the basic structure of the uncovered region M(D') — ^ 

^ Recall that V' = T> \ Q where Q is the C most expensive disks in OPT. 
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Figure 2; The general picture of the substructures in a block. The red points are the grid points 
of small squares. Dash green disks are what we have selected in %. There are five substructures in 
the block. 

For ease of notation, we use U(?^) to denote the uncovered region ]R(T>') — Figure 2 shows 

an example. Intuitively, the region consists of several “strips” along the boundary of %. Now, we 
define some notions to describe the structure of those strips. 

1. (Arcs) Consider a disk D G 2?' and suppose the center of D is in the square F. Let be the 
square gadget Gg(F), and without loss of generality assume the line DgDt is horizontal and 
divides the plane into two halfplanes and H~. D may contribute at most two uncovered 
arcs, one in H'^ and one in H~. Let us first focus on 22'^. 9D intersects dT-L at several points 
® in . The uncovered arc is the segment of 9D starting from the first intersection point 
and ending at the last intersection point. ® We can define the uncovered arc for H~ in the 
same way (if \dDr\dT-L\ ^ 0). Figure 3 illustrates why we need so many words to define an arc. 
Essentially, some portions of an arc may be covered by some other disks in 22, and the arc 
is broken into several pieces. Our definition says those pieces should be treated as a whole. 
In this paper, when we mention an arc, we mean an entire uncovered arc (w.r.t. the current 
22). Note that both endpoints of an arc lie on the boundary of 22. 

2. (Subarcs) For an arc a, we use a[A, B] to denote the closed subarc of arc a from point A 
to point B. Similarly, we only write a{A,B) to denote the corresponding open subarc (with 
endpoints A and B excluded). 

3. (Central Angle) Suppose arc a is part of cID for some disk D with center D. The central 
angle of a, denoted as A{a) is the angle whose apex (vertex) is D and both legs (sides) are 
the radii connecting D and the endpoints of a. We can show that /(a) < vr for any arc a (See 
Lemma 39 in Appendix A.l) 

® The number must be even. 

® Note that an uncovered arc may not entirely lie in the uncovered region U(2^) (some portion may be covered by 
some disks inT-i). 
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Figure 3: The figure gives an example of an arc. The blue curves are part of the boundary of %. 
The red curve is an uncovered arc. 


envelope 



baseline 


Figure 4: A substructure. The baseline b consists of the red arcs which are the part of consecutive 
boundary of &H. Qs, Qt are the endpoints of b. The black curves are are uncovered arcs. The bold 
black arcs form the envelope. The arc a -< c because A ^ C and B D. 

4. (Baseline) We use dV. to denote to be the boundary of Ti. Consider an arc a whose endpoints 
Pi, P 2 are on dTi. We say the arc a cover a point P G dP, if P lies in the segment between 
Pi and P 2 along dP. We say a point P G dP can be covered if some arc in P' covers P. A 
baseline is a consecutive maximal segment of dP that can be covered. We usually use b to 
denote a baseline. 

5. (Substructure) A substructure St(b,^) consists of a baseline b and the collection A of arcs 
which can cover some point in b. The two endpoints of each arc a € A are on b and Z(a) is 
less than tt. Note that every point of b is covered by some arc in A. Figure 4 illustrates the 
components of an substructure. 

Occasionally, we need a slightly generalized notion of substructure. For a set A of uncovered 
arcs, if they cover a consecutive segment of the boundary oiP, A also induces a substructure 
denoted as St[^]. 

Arc Order: Now we switch our attention to the order of the arcs in a substructure St(b,A). 
Suppose the baseline b starts at point Qs and ends up at point Qt- Consider any two points Pi 
and P 2 on the baseline b. If Pi is closer to Qs than P 2 along the baseline b, we say that Pi appears 
earlier than P 2 (denoted as Pi ^ P 2 ). Consider any two arcs a and c in A. The endpoints of arc a 
are A and B, and the endpoints of arc c are C and D. All of points A, B, C, D are on the baseline 
b. Without loss of generality, we assume that A ^ B, C ~< D and A ^ C. If B D, We say arc 
a appears earlier than arc c (denoted as a -< c). Otherwise, we say a and c are incomparable. See 
Figure 4 for an example. It is easy to see that -< defines a partial order. 

Adjacency: Consider two arcs a (with endpoints A ^ B) and c (with endpoints C -< D). If a -< c 
and C -< B, we say that a and b are adjacent (we can see that they must intersect exactly once), and 
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c is the adjacent successor of a. Similarly, we can define the adjacent successor of subarc a[Pi, P 2 ]. 
If c is the adjacent successor of a, meanwhile c intersects with subarc a[Pi, P 2 ], we say that c is the 
adjacent successor of subarc a[Pi,P 2 \. Among all adjacent successors of a[Pi,P 2 ]) we call the one 
whose intersection with a[Pi, P 2 ] is closest to Pi the jirst adjacent successor of a[Pi, P 2 ]. 

In order to carry out the dynamic program in Section 6, we need to properly orient the baseline 
for each substructure so that the (partial) order or the arcs is well defined. Our final solution in 
each substructure can be represented as a path (which is a segment of the boundary of the union 
of chosen disks). Our dynamic program essentially needs to determine such a path. To be precise, 
we provide a formal definition of a valid path, as follows. 

Definition 7 (A Valid Path). Consider a substructure St(b,^). Suppose the baseline b is oriented 
from Qs to Qt- A valid path is a path from Qs to Qt which consists of a sequence of subarcs 
{ai\Qs,Qi], a 2 [Qi, Q2], ■■■, CLk[Qk-i,Qt]} (baseline segments are considered as subarcs as well). 
For any ai, Oj+i is its adjacent successor (so ai -< Oj+ij. Qi is the intersection points of the arcs 
Qi and Oj+i. 

Note that the baseline from Qg to Qt is a trivial valid path (we do not consider any coverage 
requirement yet). Among all the valid paths in a substructure, there is one that is maximal in 
terms of the coverage ability, which we call the envelope of the substructure. 

Definition 8 (Envelope of a Substructure). Consider a substructure St(b,.A). Suppose the baseline 
b is oriented from Qs to Qt- The envelope of St is the valid path {ai[(5s,Qi], a2[Qi,Q2], ■ ■ ■, 
o-k[Qk-i,Qt]} where Oj+i is the first adjacent successor of at for all i & [k]. 

Coverage: Consider a substructure St(b,.4.). Consider an arc a with endpoints A and B on 
baseline b. We use b[A, B] to denote the segment of b that is covered by a. We say that the region 
surrounded by the arc a and b[A, B] is covered by arc a and use M(a) to denote the region. We note 
that the covered region M(o) is with respected to the current TL. Similarly, consider a valid path 
Path. The region covered by Path is UaePathl^(fl) (the union is over all arcs in Path) and is denoted 
by M(Path). Finally, we define the region covered by the substructure St, denoted by M(St), to be 
the region covered by the envelope of St. 

5 Simplifying the Problem 

The substructures may overlap in a variety of ways. As we mentioned in Section 2, we need to 
include in P more disks in order to make the substructures amenable to the dynamic programming 
technique. However, this step is somewhat involved and we decide to postpone it to the end of the 
paper (Section 7). Instead, we present in this section what the organization of the substructures 
and what properties we need after including more disks in P for the final dynamic program. 

Self-Intersections: In a substructure St, suppose there are two arcs a and c in ^ with endpoints 
A, B and C, D respectively. liA^B^C^D and a and c cover at least one and the same point 
in V, we say the substructure is self-intersecting. In other words, there exists at least one point 
covered by two non-adjacent arcs in a self-intersecting substructure. See Figure 9 for an example. 
Self-intersections are troublesome obstacles for the dynamic programming approach. So we will 
eliminate all self-intersections in Section 7. In the rest of the section, we assume all substructures 
are non-self-intersecting and discuss their properties. 
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Lemma 9 (Single Intersection Property). For any two arcs in a non-self-intersecting substructure, 
they have at most one intersection in the interior of the substructure. 

Proof. We prove by contradiction. Suppose a, and aj belong to the same substructure. Oj and aj 
intersect at point A and B. Since the substructure is non-self-intersecting, ai and aj lie on the same 
side of line A and B. Because the two radii of a* and aj are equivalent, the sum of central angles 
of at and Oj equals 27r. Thus at least one central angle of a* and aj is no less than vr, rendering a 
contradiction to the fact that the central angle (defined in Sec. 4) of any arc is less than vr. □ 

Base on the single intersection property, we can easily get the following property. 

Lemma 10. In a non-self-intersecting substructure, if a point is covered by two arcs a and b, a is 
adjacent to b. 

Order Consistency: There are two types of relations between substructures which affect how the 
orientations should be done. One is the overlapping relation and the other is remote-correlation. 
See Figure 2 for some examples. 

As we alluded in Section 2 , the two substructures which contain different related active regions 
of the same gadget interact with each other, i.e.. 

Definition 11 (Remote correlation). Consider two substructures St^ and Sti which are not over¬ 
lapping. They contain two different active regions of the same gadget respectively (recalling that 
one gadget may have two different active regions, one in , one in H~). We say that they are 
remotely correlated or R-correlated. See Figure 2. 

There are two possible baseline orientations for each substructure (clockwise or anticlockwise 
around the center of the arc), which gives rise to four possible ways to orient both St^ and St/. How¬ 
ever, there are only two (out of four) of them are consistent (thus we can do dynamic programming 
on them). More formally, we need the following definition: 

Definition 12 (Local Order Consistency). Consider two substructures Stu(bu,MM) and St/(b/,M/). 
There is an orientation for each substructure, such that the partial orders of the disks for both 
substructures are consistent in the following sense: It can not happen that a* -< bi in substructure 
Stu but aj -< bj in St/, where ai, aj G dDa, bi,bj G dD/, and ai,bi G Au, aj, bj G Ai. 

We show in the following simple lemma that the local order consistency can be easily achieved 
for the two substructures containing different active regions of the same gadget. 

Lemma 13 (Local Order Consistency). Consider two substructures St^ and St/ which are R- 
correlated. Each of them contains only one active region. Then the local order consistency holds 
for the two substructures St„ and St/. 

Proof. We consider two substructures St„ and St/. The arcs not in the active regions have no 
influence on the order consistency since each of them only appears in one substructure. So, we only 
need to consider the order of the arcs in two active regions. We use the same notations as those on 
the RHS of Figure 1. We orient the upper baseline from Qg to Qt, and the lower baseline from Pg 
to Pt. Suppose two arcs ci, C 3 belong to disk D^, and two arcs C 2 , C 4 belong to the disk D^. Assume 
C 4 -< C 3 in substructure St^, and ci -< C 2 in substructure St/. There must exist another intersection 
point on each side of the line connecting the two intersections of (ci,C 2 ), ( 03 , 04 ). This contradicts 
the fact that two unit disks have at most two intersections. □ 
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Then, we discuss the situation where two substructures overlap. We first need a few notations. 
For an arc a, we use D(a) to denote the disk associated with a. For a substructure St(b,^), we let 
D(b) be the set of disks that contributes an arc to the baseline b. 

Definition 14 (Overlapping Relation). Consider two substructures Sti(bi,^i) and St 2 (b 2 ,^ 2 ) and 
the point set V. We say Sti and St 2 overlap when there are two arcs a ^ Ai and b £ A 2 such that 
both a and b can cover at least one and the same point of V . 

The dynamic program need the following property. 

Proposition 15. The orientations of the two overlapping substructures should be different (i.e., if 
one is clockwise, the other should be anticlockwise). 

As different substructures may interact with each other, we need a dynamic program which can 
run over all substructures simultaneously. Hence, we need to define a globally consistent ordering 
of all arcs. 

Definition 16 (Global Order Consistency). We have global order consistency if there is a way to 
orient the baseline of each substructure, such that the partial orders of the disks for all substructures 
are consistent in the following sense: It can not happen that a* -< bi in substructure Stj(bj,^j) but 
aj -< bj in Stj(bj,^j), where ai,aj £ dDa-,bi,hj £ 9Db and ai,bi £ Ai,aj,bj £ Aj. 


Substructure Relation Graph ©: we construct an auxiliary graph ©, called the substructure 
relation graph, to capture all R-correlations and overlapping relations. Each node in & repre¬ 
sents a substructure, ff two substructures are R-correlated, we add a blue edge between the two 
substructures. If two substructures overlap, we add a red edge. 

Consider a red edge between Sti(bi,^i) and St 2 (b 2 ,.A 2 ). If baseline bi is oriented clockwise 
(around the center of any of its arc), then b 2 should be oriented counterclockwise, and vise versa. 
The blue edge represents the same orientation relation, i.e., if Sti(bi,^i) and St 2 (b 2 ,.A 2 ) are R- 
correlated, bi and b 2 should be oriented differently. 

It is unclear how to orient all baselines if G is an arbitrary graph. So we need to ensure that & 
has a nice structure. 

Definition 17 (Acyclic 2-Matching). IFe say the substructure relation graph 6 is an acyclic 2- 
matching, if G is acyclic and is composed by a blue matching and a red matching. In other words, 
G contains only paths, and the red edges and blue edges appear alternately in each path. 

If © is a acyclic 2-matching, we can easily assign each substructure a global arc order consistent 
orientations. 

Point-Order Consistency: Similarly to the arc order consistency, we also need define the point- 
order consistency, which is also crucial for our dynamic program. 

Definition 18 (Point Order Consistency). Suppose a set Vco of points is covered by both of two 
overlapping substructures Sti(bi,Mi) and St 2 (b 2 ,M 2 ). Consider any two points Pi,P 2 £ Vco and 
four arcs 01,02 £ Mi, 61,62 £ A 2 . Suppose Pi £ M(ai) n M( 6 i) and P 2 £ M(a 2 ) nM( 62 )- But 
Pi ^ M(a 2 )UM( 62 ) and P 2 ^ M(ai)UM( 6 i). We say Pi and P 2 are point-order consistent if ai 02 
in Sti and bi -< 62 in St 2 . We say the points in Vco satisfy point order consistency if all pairs of 
points in Vco are point-order consistent. 
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After introducing all relevant concepts, we can finally state the set of properties we need for 
the dynamic program. 

Lemma 19. After choosing %, we can ensure the following properties holds: 

PI. (Active Region Uniqueness) Each substructure contains at most one active region. 

P2. (Non-self-intersection) Every substructure is non-self-intersecting. 

P3. (Acyclic 2-Matching) The substructure relation graph 6 is an acyclic 2-matching, i.e., S 
consists of only paths. In each path, red edges and blue edges appear alternately. 

Pf. (Point Order Consistency) Any point is covered by at most two substructures. The points 
satisfy the point order consistency. 

How to ensure all these properties will be discussed in detail in Section 7. Now, everything is 
in place to describe the dynamic program. 

6 Dynamic Programming 

Suppose we have already constructed the set % such that Lemma 19 holds (along with an orientation 
for each substructure). Without loss of generality, we can assume that the remaining disks can cover 
all remaining points (otherwise, either the original instance is infeasible or our guess is wrong). In 
fact, our dynamic program is inspired, and somewhat similar to those in [2,20,27]. 

DP for Two Overlapping Substructures: For ease of description, we hrst handle the case 
where there are only two overlapping substructures. We will extend the DP to the general case 
shortly. Suppose the two substructures are Sti(bi,^i) and 812 ( 02 ,^ 2 ), bi is oriented from Pg to 
Pt and b 2 is oriented from Qg to Qt 

A state of the dynamic program is a pair <1> = (P, Q) where P is an intersection point of two 
arcs in substructure Sti and Q is an intersection point of two arcs in substructure St 2 . Fix the state 
‘h = {P,Q) and consider Sti. Let bp and tp be the two arcs intersecting at P. Suppose bp -< tp 
with endpoints [A, B), (C, D) respectively. We call arc bp the base-arc and tp the top-arc for point 
P. ' Our DP maintains that the base-arc is always paid in the subproblem, and intermediate state. 

Given state <I> = (P, Q), now we describe the subproblem associated with the state. Again first 
focus on Sti. Intuitively, a feasible solution to the subproblem restricted to Sti (resp. St 2 ) is a 
valid path starting from point P (resp. Q) and terminating at Pt (resp. Qt). More specifically, we 
construct a substructure St) (bi[P],Mi[P]) (See Figure 5): 

• bi[P] is the concatenation of subarc bp[P,B] and the original baseline segment bi[B,Pt]. All 
arcs in bi[P] have cost zero. 

• Mi[P] consists of all arcs a' G Ai such that bp -< a' (of course, with the portion covered by 
bi[P] subtracted). The cost each such arc is the same as its original cost. 

^ If P is the tail endpoint of an arc (so P is on the baseline), P only has a base-arc (no top-arc), which is the 
baseline arc it lies on. 
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Figure 5: The figure explains the dynamic program of two overlapping substructures. The left 
figure shows the subproblem OPT(P, Q). The goal of OPT(P, Q) is to find minimum valid paths 
for PPt and QQt respectively in set Ai [P] U A 2 [Q] such that the paths cover all points of P [P, Q]. 
The right figure illustrates one of its four smaller subproblems OPT(P*,Q). 


Similarly, we consider St 2 and the intersection point Q, and construct St2[(5] with baseline h 2 [Q] 
and arc set A 2 [Q]- We use P(a) (or 'P{A)) to denote the points covered by a (or (w.r.t. the 
original baseline). Let the point set P[P, Q] that we need to cover in the subproblem <1>(P, Q) be 


P[P, Q] = V{Ai[P]) + V{A2[Q]) - V{hp) - V{hQ). (1) 


We note that the minus term —P{bp) — P{bQ) is not vacuous as bp (resp. may cover some 
points in ^ 2 [Q] (resp. ^i[P]), and it is important that we do not have to cover those point (this 
subtlety is crucial in the correctness proof of the DP). The goal for the subproblem <h(P, Q) is to 
find two valid paths with minimum total weight, one from P to Pt and one from Q to Qt, such 
that they together cover all points in P[P, Q]. Note that the weights of both base-arcs bp and bg 
should be included in any feasible solution as well. 

Suppose bp{P,B] intersects its first successor at P^ (called base-adjacent point) and tp{P,D] 
intersects its first successor at P* (called top-adjacent point). Similarly, we dehne Q^, Q* in St 2 in 
exactly the same way. 

Now, computing the optimum for subproblem d>(P, Q) reduces to computing the optima for 
four smaller subproblems OPT(P^,Q), OPT(P*,Q), OPT(P, Q^) and OPT(P, Q*). We define two 
Boolean variables Ip (reps. Ig) indicating whether we can move from Vpb g to Vp^g without 
choosing a new arc. Formally, if Pp,Q = Vpb g, Ip = 0. Otherwise, Ip = 1. Similarly, if Vp^g = 
Vpgb, Ig = 0. If not, Ig = 1. The dynamic programming recursion is: 


OPT(P, Q) = min < 


OPJ{P,Q’^) + Ig-oo, 
OPT {P,Q^) + w[bg], 
OPT(P^Q) + /p•(X), 
OPT(P*,Q) + u;[5p], 


add no new arc ; 
add base-arc bg] 
add no new arc ; 
add base-arc bp. 


( 2 ) 


The optimal value we return is OPT(Ps,(5s). Now, we prove the correctness of the DP in the 
following theorem. We note that both the point-order consistency and the subtlety mentioned 
above play important roles in the proof. 
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Theorem 20. Suppose that we have two overlapping substruetures Sti(bi,^i) and St 2 (b 2 ,^ 2 )- 
Further suppose that bi and b 2 are oriented in a way such that the point-order consistency holds. 
Then, the cost of the optimal solution equals to OPJ{Ps, Qs) (which is computed by (2)). 

Proof. Consider subproblem OPT(P, Q). As we know the optimal solution of OPT(P, Q) should be 
two valid paths. One is from P to Pt and the other is from Q to Qt. Suppose they are Pathi = 
{ai[P,Pi], a 2 [Pi,P 2 ], •••, ai[Pi-i,Pi\, ak[Pk-i,Pt]} and Path 2 = {bi[Q,Qi], b 2 [Qi,Q 2 \, 
bi[Qi-i,Qi], bi[Qi-i,Qt]}. We can see that it suffices to prove that at least one of the two 
statements is true. 

• The pair of paths (Pathi—{ai[P, A)}, Path 2 ) is the optimal solution to subproblem OPT(Pi, Q). 

• The pair of paths (Pathi, Path 2 —{6i[(5, Qi)}) is the optimal solution to subproblem OPT (P, Qi). 

We prove by contradiction. Assume that both of the above statements are wrong. Suppose bp and 
bq are the base-arcs for state ^(P, Q), i.e., bp intersects with oi at point P and bq intersects with 
bi at point Q. We use P(Path) to denote the point set which is covered by path Path. Recall that 
P[P, Q] = P{Ai[P]) + V{A 2 [Q]) — V{bp) — V{bq). Since Pathi U Path 2 is the optimal solution for 
OPT(P, Q) (hence feasible), we have that V[P,Q] = P(Pathi) -|-P(Path 2 ) —V{bp) — V{bq). Then, 
the pair of paths (Pathi — {ai[P, Pi)}, Path 2 ) is the optimal solution for the subproblem in which 
we need to pick one path from Pi to Pt and one from Q to Qt to cover the points in 

P(Pathi) + P(Path 2 ) - V{bp) - V{bq) - P(ai) = P[P, Q] - P(ai). 

If not, we can get a contradiction by replacing Pathi — ai[P, Pi] and Path 2 with the optimal solution 
of the above subproblem, resulting in a solution with less weight than OPT(P, Q). Since the hrst 
statement is wrong, we must have that 

P[P,Q]-P(ai)/P[Pi,Q] 

(otherwise (Pathi — {ai[P, Pi)}, Path 2 ) is optimal for ^{Pi,Q)). We note that the LHS C RHS. 
Plugging the definition (1), we have that 

V{Ai[P\) + V{A 2 [Q]) - V{bp) - V{bq) - P(ai) / P(Ai[Pi]) + V{A 2 [Q\) - P(ai) - V{bq). 

A careful (elementwise) examination of the above inequality shows that it is only possible if 

P(6p)n(P(A2[Q])-P(5Q))/0. 

Repeating same argument, we can see that if the second statement is wrong, we have that 

P(6Q)n(P(Ai[P])-P(6p))/0. 

Hence, there exist bi G A.2[<5] and Oj G A.i[P] such that V{bp) ^^V{bi) ^ 0, and V{bq) ^^V{aj) A 0- 
However, this contradicts the point-order consistency because of 5p -< Oj and bq <bi. 

Thus, one of the two statements is trne. W.l.o.g, snppose (Pathi — Hi]}, Path 2 ) is 

the optimal solution to snbproblem OPT(Pi,Q). If Pi is the top-adjacent point of P, through 
OPT(P, Q) = OPT(Pi,Q) -|- w\bp\ in DP(2), we can get the optimal solution of subproblem 
OPT(P, Q). If not, suppose P^,Q^ is the top-adjacent point of P,Q, using the same argument, we 
can prove that (Pathi, Path 2 ) is the optimal solution to at least one of the subproblems OPT(P, Q*) 
or OPT(P*,Q). Thus, we can get the optimal solution for 0PT(P, Q) by our DP. 

□ 
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DP for the general problem: Then, we handle all substructures together. Our goal is find a 
valid path for each substructure such that minimizing the total weight of all the paths. We can 
see that we only need to handle each path in the substructure relation graph S separately (since 
different paths have no interaction at all). Hence, from now on, we simply assume that © is a path. 

Suppose the substructures are {Stfc(bfc,^fc)}fce[m]- We use and to denote two endpoints 
of bfc. Generalizing the previous section, a state for the general DP is = {Pk}ke.[m]-: where Pk is 
an intersection point in substructure Stfc. We use bp^,, tp^,, Pj^, P^ to denote the base-arc, top-arc, 
base-adjacent point, top-adjacent point (w.r.t. Pk) respectively. For each k G [m], we also define 
St^^'‘\bk[Pk], Ak[Pk]) in exactly the same way as in the previous section. Let P [{Pfc}fce[m]] be the 
point set we need to cover in the subproblem: 

P [{Pk}k^[n.]] = U PiMPk]) - U Pibpk). 

kG[m] kG[m] 

The subproblem OPT({Pfc}fcg[^]) is to find, for each substructure St^, a valid path from Pk to Bk, 
such that all points in P[{Pk}k&[m]\ can be covered and the total cost is minimized. 

The additional challenge for the general case is caused by R-correlations. If two arcs (in two 
different substructures) belong to the same disk, we say that they are siblings of each other. If we 
processed each substructure independently, some disks would be counted twice. In order to avoid 
double-counting, we should consider both siblings together, i.e., select them together and pay the 
disk only once in the DP. 

In order to implement the above idea, we need a few more notations. We construct an auxiliary 
bipartite graph *B. The nodes on one side are all disks in P' \'H, and the nodes on the other side 
are substructures. If disk Dj has an arc in the substructure Stj, we add an edge between Dj and 
Stj. Besides, for each arc of baselines, we add a node to represent it and add an edge between the 
node and the substructure which contains the arc. Because the weight of any arc of baselines is 
zero, it shall not induce contradiction that regard them as independent arcs. In fact, there is a 1-1 
mapping between the edges in IB and all arcs. See Figure 6 for an example. 

Fix a state <I> = {Pk}k&[m]- For any arc a in Stk (with intersection point Pk and base-arc bpk), 
a has three possible positions: 

1. a ^ bpk- we label its corresponding edge with “unprocessed”; 

2. a = bpk'- we label its corresponding edge with “processing”; 

3. Others: we label its corresponding edge with “done”. 

As mentioned before, we need to avoid the situation where one arc becomes the base-arc first (i.e., 
being added in solution and paid once), and its sibling becomes the base-arc in a later step (hence 
being paid twice). With the above labeling, we can see that all we need to do is to avoid the states 
in which one arc is “processing” and its sibling is “unprocessed”. If disk D is incident on at least 
one “processing” edge and not incident on any “unprocessed” edge, we say the D is ready. Let 
TZ be the set of ready disks. For each ready disk D, we use Ap(D) to denote the set of neighbors 
(i.e., substructures) of D connected by “processing” edges. We should consider all substructures in 
A^p(D) together. 

Again, we need in our DP indicator variables to tell us whether a certain transition is feasible: 
Formally, if P[{Pk}k^[m^ = P[\Pk][Pi]{i\\, let A = 0. Otherwise, let A = 1. For ease of notation, 
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Disks 



Figure 6: The bipartite graph which is used for marking the ready disks. The nodes on upper side 
represent the disks. The nodes on the lower side represent the substructures. If Dj has an arc in 
Stj, we add an arc between them. 


for a set {ek}k&{ra] and 5 C [m], we write [ek][e[]s = {ek}ke[m]\s U {e'jieS- Hence, 

[Pk][Pt]{i} = {Pk}ke[m]\i U Pi and 


[PkM ]Afp(D) — {Hfc}fce[m]\Arp(D) U {Pi}i£Np{D) 
Then we have the dynamic program as follows: 


OPT ({-Pfc}fcG[m]) 


= mm 


{OPT([Pfc][J’']{,|)+I,.oo}, 
minDe7^{0PT ([Hfc][T’%p(D)) +w^d} , 


add no disk 
add disk D 


(3) 


Note that in the second line, the arc(s) in iVp(D) are base-arcs (w.r.t. state 

In the rest of the section, we prove the correctness of the dynamic program. If we use the solution 
of a smaller subproblem OPT(<I>') to compute subproblem OPT($), we say <I> can be reached from 
(denoted as <!>' —)> ^h). If <I> can be reached from initial state <I>o({Hfc | Pk = Bk}k&[m])j we say 
the state is reachable, which is denoted by ^>o —?■ 

We start with a simple consequence of our DP: there is no double-counting. 

Lemma 21. If two arcs in the solution belong to the same disk, their weights are counted only once 
in the DP (3). 

Proof. From the DP, we can see that the weight of an arc is counted only when it becomes a base- 
arc (its label changes from “unprocessed” to “processing”). If the two sibling arcs a, b (belonging to 
disk D) in the solution are counted twice, there exist two states and $2 such that (l)d>i —?■ $ 2 ; 
(2)a is a base-arc in <I>i, but b is not a base-arc in <l>i, (3) 5 is a base-arc in <I' 2 . So, in <I>i or any 
state before that, arc b is “unprocessed”. However, a can become a base-arc only when D is ready, 
rendering a contradiction. □ 

Now, we prove the correctness of the dynamic program. The proof is a generalization of Theo¬ 
rem 20. 
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Theorem 22. Suppose & is a path. All baselines are oriented such that all properties in Lemma 19 
hold. Then, the optimal cost for the problem equals to OPT({Tfc};jg[m]) (computed by our DP (3)j. 

Proof. Suppose the set {Pathfcj^jgj^] of paths is the optimal solution. We need to prove (1) the 
hnal state d>({Tfc}fcg[^]) is reachable, (2) for any reachable state = {Pk}ke[m]y OPT({Pfc};i,g[^]) 
computes the optimal solution for the corresponding subproblem (that is to find one valid path 
from Pk to Bk for each substructure St^ to cover all point in V[{Pk}k£[m]]j such that the total 
cost is minimized). 

We first prove the hrst statement. Suppose the state d> is reachable, we prove it can reach 
another state if <1> is not the final state (i.e., we do not get stuck at <h). Let us prove it by 
contradiction. Assume we get stuck at state That means there is no ready disk in $. Note 
that each substructure, say St, is incident on exactly one “processing” arc, say arc a (which is the 
base-arc in St), a’s sibling, say b (in SL), must be labeled “unprocessed” (otherwise the disk would 
be ready). Consider the base-arc (or “processing” arc), say a', in Sth So we have a' ^ 6 in Sth 
Again the sibling b' of a' must be an “unprocessed” arc in St. ® So we have 6 ' -< a in St, which 
contradicts the global arc-order consistency. 

Now, we prove the second statement. We consider state $ = {Pk}k&[m]- Suppose the opti¬ 
mal solution for subproblem ^{Pk}ke[m] is the set of paths PS = {Pathfej^gj^j where Pathfc = 
(ofc^, 0^2 ■ • •, afc„)- Consider the states <hArp(D) := [-PA:][-F’i*]iVp(D) for all D gTZ. Obviously, we can see 
from our DP that <hArp(D) Define for each D G 7^, a set of paths 

PS 7 Vp(D) = {Pathfc}fcg[m]_ 7 Vp(D) U {Pathj - {aji}}jgjVp(D)- 

It suffices to prove that there exists at least one D G TZ snch that PSjv (d) is the optimal solution 
for OPT($jVp(D))- 

Consider a substructure St*. Suppose the intersection point in St* of is Pi and the base-arc at 
point Pi isbpi. For each i G [m], let St 7 r(i) be the only substructure (if any) overlapping with St*. So 
is the base-arc of St.n-(j). Using exactly the same exchange argument in Theorem 20, we can 
show that if PS 7 Vp(D) is not the optimal solution for 0 PT(<l> 7 Vp(D))) there exists some i G 7Vp(D) such 
that Si happens, where Si is the following event: there exists an arc in St 7 r(i) with bp^^i^i-^ 

snch that 

iP(6p,)nP(/3^(i))/0. 

We use Si to denote the above event. Thus if there is no D G 7^ such that PS 7 Vp(D) is the optimal 
solution for 0 PT(<I>^p(D)), we have 

/\ [ V f.j=True. (4) 

Np{D)&n \ieAfp(D) / 

Converting the conjunctive normal form (CNF) to the disjunctive normal form (DNF), we get 


V 


f’fc. I = True, 


(fci,...,fc|.jj^|)gnog7jAfp(D) yis|7?.| 


®To see that b' is in St, note that S is a path and St is only R-correlated with St' (and vice versa). 
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where nDe'R,A^p(D) means the Cartesian product of all iVp(D) in TZ. We call each A^e| 7 ^| a clause 
(note that ki indexes a substructure). If we can prove that every clause is false, then obviously, (4) 
is wrong, resnlting in a contradiction. 

Now, we show that every clause is false. First, we consider the case that both end nodes of S 
are incident to red edges. W.l.o.g., suppose the two nodes of S are Sti and St 2 . Thus, they are 
R-correlated with other substructures. We know if one substructure Stj is not R-correlated with 
others, every clause must contain the corresponding event A (since the disk corresponding to its 
base-arc must be ready and in TV). Hence, every danse contains £i and £ 2 - Moreover, for each 
pair (Stj, Stj) of R-correlated snbstructures, each clause should contain either A or £j. Suppose 
the length of the path © is Because red and blue edges alternates, there are R-correlated 
snbstructnre pairs and overlapping substructure pairs. We should select -|- 2 terms in each 
clause. Because + 2 > there exists a pair of overlapping substructures (Stj,St^(j)) such 
that both £i and appear in the clause. To make the clause true, we must have 

£i = CPibpi) n / 0) = True and n (5i ^ = True, 

where bp^^^i) and j3i)~ bpi- It yields a contradiction to the point-order consistency. 

Next, we consider the remaining case where at least one end of the path S is a blue edge, 
meaning the substructure on the end does not overlap with any other substructure. W.l.o.g., 
suppose the end node is Sti and it is R-correlated with St 2 . The event £i = Vihpi) n'P(/3^(i)) ^ 0) 
is always false since /3jr(i) does not exist. So all clause containing £i is false. To make each of 
the remaining clauses true, £2 must be true, and the case reduces to the previous case (by simply 
omitting node Sti). So the same argnment again renders a contradiction. This completes the proof 
of the theorem. 

□ 


7 Constructing 71 

In this section, we describe how to construct the set TL in details. We first include in Ti all square 
gadgets. The boundary of Ti consists of several closed curves, as shown in Fignre 2. Ti and all 
nncovered arcs define a set of substructures. 

First, we note that there may exist a closed curve that all points on the curve are covered by 
some arcs (or informally, we have a cyclic substructure, with the baseline being a cycle). We need 
to break all such baseline cycles by including a constant number of extra arcs into Ti. This is easy 
after we introduce the label-cut operation in Section 7.1, and we will spell out all details then. 
Note that we cannot choose some arbitrary envelope cycle since it may ruin some good properties 
we want to maintain. 

From now on, we assume that all baselines are simple paths. Now, each closed curve contains 
one or more baselines. So, we have an initial set of well defined substructures. The main purpose 
of this section is to cut these initial substructures such that Lemma 19 holds. 

We will execute a series of operations for constructing Ti. We first provide below a high level 
sketch of our algorithm, and outline how the substructures and the snbstructnre relation graph S 
evolve along with the operations. 

• (Section 7.1) First, we deal with active regions. Sometimes, two active region may overlap 
significantly and become inseparable (formally defined later), they essentially need to be dealt 
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as a single active region. In this case, we merge the two active regions together (we do not 
need to do anything, but just to pretend that there is only one active region). We can also 
show that one active region can be merged with at most one other active region. For the 
rest of cases, two overlapping active region are separable, and we can cut them into at most 
two non-overlapping active regions, by adding a small number of extra disks in %. After 
the merging and cutting operations, each substructure contains at most one active region. 
Hence, the substructures satisfy the property (PI) in Lemma 19. Moreover, we show that if 
any substructure contains an active region, the substructure is limited in a small region. 

• (Section 7.2) We ensure that each substructure is non-self-intersecting by a simple greedy 
algorithm. After this step, (P2) is satisfied. 

• (Section 7.3) In this step, we ensure that substructure relation graph S is a acyclic 2-matching 
(P3). The step has three stages. First, we prove that the set of blue edges forms a matching. 
Second, we give an algorithm for cutting the substructures which overlap with two or more 
other substructures. After the cut, each substructure overlaps with no more than one other 
substructure. So after the first two stages, we can see that S is composed of a blue matching 
and a red matching. At last, we prove that the blue edges and red edges cannot form a cycle, 
establishing © is acyclic. 

• (Section 7.4) The goal of this step is to ensure the point-order consistency (P4). We first show 
there does not exist a point covered by more than two substructures, when © is an acyclic 
2-matching. Hence, we only need to handle the case of two overlapping substructures. We 
show it is enough to break all cycles in a certain planar directed graph. Again, we can add a 
few more disks to cut all such cycles. 

• (Section 7.5) Lastly, we show that the number of disks added in T-L in the above four steps is 
0{K‘^), where K = ^ and L and /r are side lengths of block and small square respectively. 

7.1 Merging and Cutting Active Regions 

If two active regions overlap in the same substructure, we need to either merge them into a new 
one or cut them into two non-overlapping new ones. As we know, each gadget may have two active 
regions. Suppose active regions Ari and Ar 2 belong to the same gadget Gg, while Ar)^ and Ar 2 
belong to a different gadget Ggk Due to R-correlations, we need consider the four active regions 
together. 

First, we consider the case where Ari overlaps with Ad^, and Ar2 overlaps with Ar^. We need 
the following important concept order-separability, which characterizes how the two sets of arcs 
overlap. 

Definition 23 (Order-separability). Consider a substructure St(b,A). Ai,A 2 are two disjoint 
subsets of A. If Ai,A 2 satisfy that 

a <b, for any a G Ai and b G A 2 , (5) 

we say that Ai,A 2 are order-separable. 

We use Ai , A!i, A 2 , A2 to denote the set of arc associated with active regions Ari, Ar'^, Ar2, 
Ar2, respectively. If Ai and A!i are not order-separable, we say the pair (Ari, Ar'^) is a mixture. If 
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both of (Ari,Ar'^) and (Ar 2 ,Ar 2 ) are mixtures, we say the two pairs form a double-mixture. When 
they are double-mixture, we merge them simultaneously. It only means that we regard the two 
active regions (Ari, Ar^^) as a new single active region, and (Ar 2 , Ar^) as another single active region. 

To show an active region cannot grow unbounded, we prove that the merge operations do 
not generate chain reactions. The rough idea is that if two active regions form a mixture, their 
corresponding small squares must be adjacent and their core-central areas must overlap. Due to 
the special shape of core-central areas (a narrow spindle shape), the overlapping can only happen 
between two of them, not more. 

Lemma 24. Consider two non-empty small squares T, T'. Suppose the square gadgets in the 
two squares are Gg(r) = (D^jD^) and Gg'(r') = (D'„Dj ). The active region pairs (Ari, Ar 2 ) and 
(Ar;,Ar' ) are associated with gadget Gg(r) and Gg'(r') respectively. (Ari, Ar 2 ) and (Ar'^,Ar 2 ) form 
a double-mixture. Then, the following statements hold: 

1. Their corresponding squares T, T' are adjacent; 

2. The core-central areas of Gg{T) and Gg'(r') overlap; 

3. The angle between DgDt and D'gD[ is 0(e) 

j. None of the two core-central areas can overlap with any small squares other than T and T'. 

Next, we consider the case where only one of (Ari,Ar^^) and (Ar 2 , Ar^) is a mixture. However, 
we show that it is impossible as follows. We use notation Ar(^) to denote an active region with 
arc set A. 

Lemma 25. Suppose active region pairs (Ari(.4i), Ar2(.4.2)) and (Ar)^(.4']^), Ar2(.42)) are associated 
with gadget Gg and Gg' respectively. If Ai and A!i are order-separable, then A 2 and AI 2 are also 
order-separable. 

The proofs of Lemma 24 and Lemma 25 are elementary planar geometry and we defer them to 
Appendix A.2. 

Cutting Overlapping Active Regions: After the merging stage, even if any two active regions 
overlap in the same substructure, their arcs are order-separable. We define label-cut operation as 
below to further separate them such that the baselines of all the active regions are non-overlapping. 

We consider two overlapping active regions Ari(Ai) and Ar 2 (A 2 ) in a substructure St(A). Sup¬ 
pose Ai,A 2 C a and Ai,A 2 are order-separable. We can cut St into two substructures Sti and 
St 2 with disjoint baselines, such that M(Ari) C M(Sti) and M(Ar 2 ) C M(St 2 ). In other words, if 
we assign arcs in Ai one kind of label and arcs in A 2 a different kind of label, after the label-cut 
operation, the arcs with different labels belong to different new substructures. We define label-cut 
formally as follows. 

Definition 26 (Label-cut). Consider a substructure St(b, A). Suppose A)^,A 2 are two subsets of 
A. A!i and AI 2 are order-separable. We can add two consecutive arcs in the envelope of A into 
TL. After that, St can be separated into two new substructures Sti(bi, Ai), St 2 (b 2 , A 2 ), such that 
subarcs in A'^ only belong to Ai, and subarcs in A 2 only belong to A 2 . 

Lemma 27. Consider a substructure St(b, A) and two subsets A[, A 2 of A. There exists a label-cut 
when A!i and A 2 are order-separable, and we can perform the label-cut in polynomial time. 
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envelope 



Figure 7: The example of label-cut. The left hand side illustrates the whole substructure before 
cutting. The arcs have two different labels. One is green and the other is brown. The bold black 
subarcs are what we select in the envelope. The the right hand side illustrates that each of the two 
separable substructures induced by the label-cut operation only contains arcs with the same label. 



Figure 8: The farthest disk pair of square T is (D^, D^). Suppose and intersect at point P, Q. 
The disk D(P, 2 ) is tangent to D* and at point Qs and Qt respectively. D(Qs,l) and D(Qt,l) 
intersect at D. The active region of F in is totally covered by Dom(r'’'). 


Figure 7 illustrates the process of construction. We defer its proof in Appendix A.2. 

After the label-cut operation, in each substructure, the baselines for all active regions are 
not overlapping. Thus, if any substructure contains more than one active region, consider any 
two of them, say Ai and A 2 - We can add into % one arc a along the envelope which satisfies 
ai -< a -< a2,Voi G ^1,02 £ A2- After the addition of arc a, and A2 are separated into two 
different new substructures. Repeat the above step whenever one substructure contains more than 
one active region. This establishes the active region uniqueness property (PI). 

Limiting the size of substructure which contains an active region: Now, we discuss how 
to make substructure which contains an active region bounded inside a small region. This property 
is particular useful later when we show the substructure relation graph © is acyclic. 

Suppose the gadget of square F is (D*, Dt). The line DgDt divides the plane into two halfplanes 
and H~ . and intersect at point P and Q. The boundary of disk D(P, 2) is tangent 
to Ds and at point Qs and Qt respectively. D((55,l) and D(Qt,l) intersect at point D. See 
Figure 8. We call D the dome-point of gadget Gg(r) and use Dom(r+) to denote the region 
{D{DA))-Ds-Dt)nH+. 

Lemma 28. Consider the substructure St which contains an active region of square F. The sub¬ 
structure can be cut into at most three smaller substructures, by doing label-cut twice. At most 
one of them contains the active region. Moreover, this new substructure (if any) is bounded by the 
region Dom(r’^) (or Dom(F“)j associated with F. 

We defer the proof to Appendix A.2. 
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There may exist a substructure which contains a merged active region (i.e. an active region 
which is the union of two initial active regions). Based on Lemma 24, the arcs of the two initial 
active regions belong to two adjacent small squares. Suppose the two squares are T, T', with square 
gadgets (Ds,Dt) and (0^,0^) respectively. The dome-points of T and T' in are respectively 
D and D'. We apply the operations in Lemma 28 for each initial active region. Since the angle 
between DgDt and D'gD[ is 0(e), obviously, the substructure containing the merged active region 
is small as well, i.e. bounded in Dom(r’^) U Dom(r'+). 

Since each substructure which contains an active region is small enough, we can show the 
following lemma, which will be useful for proving the acyclicity of the substructure relation graph 
6 . 

Lemma 29 (Highly parallel arcs). Suppose substructure St(b,.4.) contains an active region. The 
central angle of any arc in St is no more than 0(e). Meanwhile, there exists a line I such that the 
angle between I and the tangent line at any point of any arc a ^ A is at most 0(e). 

Actually, for the initial active region, the line DgDt satisfy the property. For the merged active 
region, we know the angle between DgDt and D'gD[ is 0(e) based on Lemma 24. Hence, we can 
still see that the line DgDt satisfies the property. 

To summarize, after all operations in this subsection, we can ensure that each substructure con¬ 
tains at most one active region (i.e., (PI) in Lemma 19). Moreover, we have that each substructure 
which contains an active region is small enough (so that Lemma 29 holds). 

Handling cyclic substructnres: At the end of this subsection, we deal with the problem we 
left in the very beginning of Section 7, to break all cyclic baselines. Note that this step should be 
done in the beginning. First, we consider that case that there exists a point on the baseline which 
cannot be covered by any arc of any active region. We can include any envelope arc that covers the 
point into T-L, which is enough to break the cycle. This is essentially a label-cut and do not separate 
any single connected active region into disconnected pieces. Then we consider the case where every 
point on the baseline is covered by some arc of active regions. Note that the merge operation only 
depends on the local property of two active regions, thus, we can merge active regions even in 
“cyclic substructure”. Assume that we have done all merge operations. Based on Lemma 29, we 
know one active region is very small and thus cannot cover all points on a closed curve. We pick 
one active region. Using the same operation as Lemma 28 (do two label-cuts), we can essentially 
isolate the active region and cut the original cyclic baseline to two new baselines. 

7.2 Eliminating Self-intersections 

The goal of this part is to add a few more disks into TL so that any substructure is non-self- 
intersecting (Recall the definition in Lemma 9). Note that the substructures which contain active 
region are non-self-intersecting after the process in Section 7.1. We just need to process the sub¬ 
structures without active regions in this part. 

We use a simple greedy approach. We consider one substructure St(b, ^). Suppose the endpoints 
of b are Qs and Qt, and the envelope is Qi],. ■., ak[Qk-i, Qt]}, where Qi is the intersection 

point of Oi-i and a*. We denote the endpoints of ai on the baseline b by Pi and P/ (Note that 
Pi = Qs)- Let arc set Ai be { 01 , 02 ,. ■. ,Oi}. Then we have a set {St4-4.j]}jg[fc] of substructures, 
where Sti[-4.i] is the substructure induced by arc set Ai- We consider the arcs lying on the envelope 
one by one and check whether we should add it into TL or not. Concretely, we add D(oi+i) in Ti if 
the following condition holds: 
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(a) 


(*) 


Figure 9: The process to avoid self-intersection. The left hand side is a self-intersection substructure. 
We search from point Qs along the envelope. Let arc set Ai be {ai, 02 , ■ ■ •, Oi}. Then we have a set 
{Sti[^i]}jg[A;] of substructures. If St* is non-self-intersecting but Stj+i is self-intersecting. We add 
arc Oj+i in Ti. The right hand side illustrates the two new substructures after the cut. 


• Stj is non-self-intersecting, but Stj+i is self-intersecting, 

The addition of D(ai+i) cuts the substructure into two smaller substructures. One of which is 
certainly non-self-intersecting. The other is induced by arc set A = {ai+ 2 j ■ • • See Figure 9. 

We repeat the above process until there is no self-intersection in all substructures. Furthermore, 
we can easily prove the following nice property. The proof can be found in Appendix A.3. 

Lemma 30. In each of the above iterations, one substructure St(b,.A) is cut into at most two new 
substructures. Any original arc in A cannot he cut into two pieces, each of which belongs to a 
different new substructure. 

To summarize, we have obtained the non-self-intersection property ((P2) in Lemma 19). 

7.3 Ensuring that © is an Acyclic 2-Matching 

We discuss how to add some extra disks in H to make S an acyclic 2-matching ((P3) in Lemma 19). 
Blue edges: First we show that the set of blue edges form a matching. 

Lemma 31. Two blue edges cannot be incident to the same node. 

Proof. Before the merge operation, the set of active region pairs forms a matching. To see this, 
note that our merge operations always apply to a double-mixture (which corresponds to merging 
two blue edges into one). Moreover, any cut operation cannot break one active region into two, 
thus has no effect on any blue edge. Hence, after all merge and cut operations, the set of active 
region pairs is still a matching. □ 

Red edges: Then, we prove that any node which has more than one incident red edges can be cut 
such that each new node(i.e., substructure) only has at most one incident red edge. 

First, we prove a simple yet useful geometric lemma stating that a point cannot be covered by 
three or more substructures. Note that from now on, all substructures have no self-intersections. 

Lemma 32. We are given a substructure St(b,A). a G A and two arcs 61,62 ^ A. 1 / 61,62 cover 
the same point on a, 61,62 should belong to the same substructure. 

Intuitively, if the two disks corresponding to 61 and 62 cover the same point, they should be close 
enough such that their corresponding square gadgets overlap (which implies 61 and 62 share the 
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Figure 10: The arcs in substructures Sti and St 2 cut the envelope of St into 7 segments. The 
sequence of the labels for those segments is 0101020. The compressed label sequence is 01020. So 
we have 5 1-segments. 


same baseline). First we prove that the minimum distance between any two disks in two different 
substructure should not be too small, i.e., their overlapping region cannot be too large. The proof 
can be found in Appendix A.4. 

In fact, essentially the same proof can be used to prove that any two different substructures 
cannot both intersect with subarc whose central angle is 0(e), as in the following corollary. 

Corollary 33. Suppose St(b,^) is a substructure without any self-intersection. Consider an arc 
a £ A and two arcs 61,62 ^ A. Suppose a' is a subarc of a whose central angle is 0{e). If both 
61,62 cover some part of a', 61,62 should belong to the same substructure. 

Combining with Lemma 29, we can easily see the following lemma. 

Lemma 34. Any substructure which contains an active region cannot overlap with two or more 
different substructures. 

Then we show how to cut the substructure St(b,^) which overlaps with more than one other 
substructures. Note that such substructure does not contain an active region based on Lemma 34. 
Suppose the envelope of St is {ai[Qs,Qi],... ,ak[Qk-i,Qt]}- St overlaps with k substructures 
Stj(b2, Aif i — 1,2,. 

If Stj overlaps with St, there exists an arc a £ Stj intersecting some envelope arc of St. Thus, 
the envelope St is subdivided into several segments by those intersection points. We can label those 
segments as follows: 

• If the segment is covered by some arc in Stj, we label it as ‘i’. 

• If there is no arc in any Stj covering the segment, we label it as ‘O’. 

See Figure 10 for an example. According to the Lemma 32, we know there is no point on the 
envelope covered by two substructures. Thus the above labeling scheme is well defined. 

Traversing those segments along the envelope, we obtain a label sequence. First, for each label 
i, we identify those maximal consecutive subsequence, which consists of only letter 0 and i, and 
starts with and end with i, and replace the subsequence by a single letter i. We obtain a compressed 
sequence. In fact, each letter, say i (f / 0), in the compressed sequence corresponds to one or more 

segments labeled with either i or 0, and the first and last one must be labeled with i. We call 

the concatenation of those segments an l(ong)-segment. Of course, a letter 0 in the compressed 
sequence corresponding to a segment with label 0. Actually, the sequence is Davenport-Schinzel 
sequence of order 2 [37], because two different labels cannot intersect (because the baselines of two 


24 



substructures cannot intersect). For example, the pattern “1212” should never appear. Thus, the 
length of the compressed sequence, i.e., the number of 1 -segments, is at most 0{k). 

Now, we discuss how to cut St into several new ones based on 1-segments. Keep in mind that our 
goal is to make sure each new substructure only overlap with one substructure of {Stjjjgjfc]. The 
cut operation is again a simple greedy procedure. Consider two consecutive 1-segments. Suppose 
they are Qj-|_i], ... ,(ij[Qj—i,Qj]} and Qj+i]; •••) Ojk[Qk—iiQk]}- 

We add into Ti the last arc aj of former 1-segment and the first arc aj+i of the later 1-segment. St 
is thus cut into two new ones. Repeat the above step for all two consecutive 1-segments in order. 

We still need to show that after the cut, every new substructure overlap at most one of {Stj}jg[fc]. 
Consider one new substructure. Notice such an original arc in A can only belong to one new 
substructure. By our cut operation, all envelope arcs of the new substructure can intersect at most 
one of {Stj}jg[fc]. Hence, the new substructure can overlap at most one Stj. 

Blue edges and red edges: After the above operations, the set of all blue edges forms a matching, 
while the set of all red edges also forms a matching. To show & is an acyclic 2-matching, it suffices 
to prove & contains no cycle. So, the rest of the section is devoted to prove the following lemma. 

Lemma 35. Suppose the side length of square is p,, where p = 0(e) and the block contains K x K 
small squares, where K = ^ and Co is an appropriate eonstant. Then, after all operations stated 
in this seetion, there is no eyele in 6. 

If there is a cycle Cyc in ©, the red edges and blue edges alternate in Cyc, which correspond to 
a sequence of substructures, each containing an active region (since it is R-correlated with another 
active region). Now, we provide a high level explanation why Lemma 35 should hold. Each active 
region is associated with a small square. A small square is very small (comparing to a unit disk 
or the whole block), and an active region is also very small. We pick a point in each small square 
and each substructure. If two substructures in Cyc overlap, the two points in them are also very 
close (i.e., 0(e)). So we can pretend the two points as one point (or we just pick a point in their 
overlapping region). For each active region Ar, we connect the point in Ar and the point in the small 
square associated with Ar. Since Both the small square and the active region are very small, the 
distance between two point is about 1 — 0(e). Thus, a cycle Cyc would present itself geometrically 
as a polygon. We can show the angle between two adjacent edges of the polygon is close to vr. So 
the size of the polygon cannot be not small (it takes a lot of edges to wrap a loop). However, the 
polygon cannot be much larger than the block. By contradiction, we prove that there is no cycle 
in S. 

Now, we formally prove Lemma 35. We first prove a geometric lemma which will be useful for 
bounding the angle between the two adjacent edges of the aforementioned polygon. 

Lemma 36. Consider two substructures Sti(bi, Ai) and St 2 (bi, A 2 ) in Cyc. Suppose Sti and St 2 
overlap (there is a red edge between them). For any two arcs a G and b G A 2 , suppose that D(a) 
and D(6) overlap and their intersection points are A and B. The central angles of a,b are 9a, Ob 
respectively. The eenters o/D(a), D(5) are Da,Db. Then ZADaB (or ZADbB) is at most {9a + 0b). 

Proof. We distinguish a few cases depending on whether the intersection points A and B lie on a 
or b or none of them. All cases are depicted in Figure 11. 

• Both A and B lie on one arc (see Figure (a)(b)). W.l.o.g., suppose they lie on arc a. Obviously, 
ZADaB is no more than 6a (or 6b) . 
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Figure 11; The four kinds of cases for two overlapping arcs. 


• If one intersection is on neither a nor h (see Figure (c)), we prove that the case cannot happen. 
Suppose A is on neither of a and b. The endpoint Ai of a is covered by disk D{b) and the 
endpoint A 2 of a is covered by disk D(a). Thus the baseline bi must intersect with b2 which 
contradicts the fact that Sti and St2 are two different substructures. 

• If one intersection point is on a but not on b and the other intersection point is on b but not 
on a (see Figure (d)), it is easy to see that Bi is covered by arc b, and ^42 is covered by a 
since their baselines do not intersect. Thus, the length of arc AB (w.r.t. D(a)) is no more 
than the sum of lengths of a and b. So ZADaB is at most 9a + Oh- 

The above cases are exhaustive, thus our proof is completed. □ 

Based on Lemma 29, Lemma 36, we can prove Lemma 35 below: 

Proof of Lemma 35: Consider two substructures Sti and St2. Suppose Sti contains the active 
region Ari , and St 2 contains the active region Ar 2 . The centers of arcs of Ari and Ar 2 locate in 
small squares Fi and F 2 respectively. The square gadgets of the two squares are Gg(Ds^, and 
Gg(Ds 2 , 0 ^ 2 ) respectively. If Sti and St2 overlap, their must exist an arc a in Sti and an arc b in 
St 2 such that a and b intersect. Suppose D(o) and D(6) intersect at points A and B. The center of 
D(a) and D(6) are Da and Dh- 

Based on Lemma 29, we know both central angle of a and b are no more than 0(e). According 
to Lemma 36, the central angle of the arc AB is at most 0(e). It means angle between the tangent 
lines of D(a) and D(6) at point A is no more than 0(e). Thus, ZDaADh is at least vr — 0(e). 

We know that all disks in the same active region are centered in one small square or two 
adjacent small squares. Moreover Lemma 29 implies all disk centers should lies in one or two 
squares. Hence, the distance between (any point in) the square and (any point in) its active region 
is at least 1 — 0(e). Construct the aforementioned polygon. Consider two adjacent edges XY and 
YZ in the polygon. We consider two cases: 

1 . Y is in the intersection of two substructures Sti and St2. We can easily see that ( 1 ) \YA\ = 
0(e); (2) \XDa\ = 0(e); (3) \ZDh\ = 0(e). Hence, we can see ZXYZ is at least tt — 0(e). 
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Figure 12: Substructure Sti(bi,^i) and St2(b2,.42) are overlapping. bi starts from Qg and ends 
up with Qt and b2 starts from Pg and ends up with Pt. There are two paths forming a cycle. 

2 . y is in a small sqnare T and X and Z are in the the two substructures associated with F. 
Since both substrnctures are bound in an 0 (e) size region (by Lemma 29 ), we can see that 
ZXYZ is at least tt — 0 (e) as well. 

Hence, we can see the polygon contains at least H( 27 r/e) nodes and the diameter of the polygon is 
at least H(l/e). But this cannot be larger than the diameter of a block, rendering a contradiction. 

□ 

To summarize, we have ensured that & is an acyclic 2 -matching ((P 3 ) in Lemma 19 ). 

7.4 Ensuring Point Order Consistence 

We have ensured that the set of red edges is a matching in S. Hence, one substructure can overlap 
with at most one other substructnre. Therefore, if we can gnarantee that the points which are 
covered by any pair of overlapping substructure are order consistent, then all points in V are order 
consistence (after all, point-order consistency is defined over a pair of snbstructnres). 

Consider two overlapping substructures Sti(bi,^i) and St2(b2,.42) and a set Vco of points 
covered by .4.1 U 4.2. Suppose Pi,P2 G Vco and 01,02 G 4.1,61,62 G 42. Recall point-order 
consistency requires that the following conditions are met 

• Pi G M(oi) n M(6i) and P2 G M(a2) n 1^(62) 

• Pi ^ M(o 2) U M(62) and P2 ^ M(ai) U M(6i). 

• (oi -< 02) (61 -< 62). 

It is helpful to consider the following directed planar graph D induced by all arcs in 4 i U 42 
in the uncovered region U(P). Regard each intersection point of arcs in 4 i U 42 as a node. Each 
snbarc is an directed edge with a direction consistent with its baseline. We use 4 —)• P to denote 
that there is a directed edge between nodes A and B in D. If there is no directed cycle in D, 
we can verify that all conditions listed above hold. Indeed, suppose the condition is not satisfied, 
which means oi -< 02, 61 62, oi, 02 G 4 i and 61,62 G 42- Suppose oi, 61 intersect at ( 4 i, Pi) and 

4 i —>■ Pi, meanwhile 02,62 intersect at (42, P2) and 42 —)• P2. Since oi -< 02, there exists a path 


27 




in Sti which goes from Bi to ^2- Similarly, there exists a path in St2 from B2 to Ai. Thus, the 
two paths and ai, 02 form a directed cycle. 

So all we have to do is to break all cycles in X>. When Ti contains a cycle, we can cut the cycle 
through adding an arc on the envelope into Ti. See Figure 12 for an example. Only arcs .Ai (or in 
A2) cannot form a cycle. So if there is a cycle, the cycle must pass through the envelope of Ai and 
A.2. Moreover, based on Lemma 29 , if Sti and St2 form a cycle, either Sti or St2 does not contains 
any active region. W.l.o.g., suppose it is St2. Based on these observations, we have our algorithms 
as follows: 

Suppose the envelope of Sti is Pathi = {ai, 02,..., Ofc}. Oj and aj are the first and last arcs 
respectively which intersect St2. Suppose the arc G A2 overlaps with a*, (if there is more than 
one such arc, we select a minimal one, w.r.t. the arc ordering) and the arc bj/ G A2 overlaps with 
Uj (if there is more than one such arc, we select a maximal one). Since a* -< Uj and they do not 
satisfy point-order consistency, we have 6^/ bj' . We can see that all arcs between 6j/ and bj' cannot 
intersect with Pathi. So we can select one arc between 6j/ and bj' to add in Ti for cutting St2 into 
two. After the cut, any cycle in D can be broken. 

After cutting St2, Sti overlaps with both of the new substructures obtained from St2. Then, we 
encounter the same situation as in Section 7.3 (a node in & has two incident red edges). We can 
apply the operation in Section 7.3 to cut Sti such that the set of red edges in © is still a matching. 

7.5 The number of disks in H 

Finally, we count collectively the total number of disks that we have added in Ti. First, we add 
the square gadget for each nonempty small square in B. The number of the disks is 0 {K‘^), where 
K = Ljfj, = 0 (l/e^). In order to cut overlapping active regions, the number of disks we add in B 
is bounded by the number of active regions. Since there are O(iF^) gadgets, we add O(iF^) disks 
in Section 7 . 1 . In Section 7 . 2 , to ensure that each substructure is non-self-intersecting and contains 
at most one active region, we design a greedy algorithm. We can see that each arc we added in B 
covers at least one intersection point of two disks in B', where B' is the set B before this step. The 
algorithm guarantees that each arc which we add in B do not cover the same intersection point on 
the boundary of B'. Since the union complexity of unit disks is linear [ 32 ] and B' contains at most 
0 {K'^) disks, there are at most 0 {K'^) intersection points on dB'. So, we add at most 0 {K^) in 
this step. In Section 7 . 3 , we break the cycles in S. The number of disks we add B is proportional 
to the number of substructures. So, again, we add at most 0 {K‘^) disks. Similarly, in Section 7 . 4 , 
we also add at most 0 {K^) disks. To summarize, we have added at most 0 {K^) disks in B. 


8 Time Complexity 


The time complexity contains three parts. The first part is to enumerate all combination of Q. 
We set C = O(^) (since we need C > \B\/e). Since K = O(^), the number of combinations 
is bounded by inP = The second part is the construction of set B. It is easy to see 

that the time cost for each operation (i.e., label-cut) is no more than O(n^). Thus, the time 
cost is 0 {K^‘n?) = 0 (n^/e^). The last part is the dynamic program. There are at most O(iF^) 
substructures and at most O(n^) intersection points in each substructure. Thus, the number of 


total states is at most 0((n 


2\K^ 


For each recursion, the time cost is at most 0 {n). Thus, the 


overall time complexity of the dynamic program is 0 {n^^ ^^) = \ 
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Overall, the total time cost is • max{n^/e'^, This finishes the proof 

of Lemma 6 . 

9 Applications 

The weighted dominating set problem (MWDS) in unit disk graphs has numerous applications in 
the areas of wireless sensor networks [17]. In this section, we show that our PTAS for WUDC can 
be used to obtain PTASs for two important problems in this domain. 

9.1 Connected Dominating Set in UDG 

The goal for the minimum-weighted connected dominating set problem (MWCDS) is to hnd a 
dominating set which induces a connected subgraph and has the minimum total weight. Clark 
et al. [11] proved that MWCDS in unit disk graphs is NP-hard. Ambiihl et al. [ 2 ] obtained the 
first constant factor approximation algorithm for MWCDS(the constant is 94). The ratio was 
subsequently improved in a series of papers [13,20,30]. The best ratio known is 7.105 [17, pp.78]. 

One way to compute an approximation solution for MWCDS is to first compute minimum 
weighted dominating set (MWDS) and then connect the dominating set using a node-weighted 
Steiner tree (NWST) [30,41]. The optimal MWDS value is no more than the optimal MWCDS 
value. After zeroing out the weight of all terminals, the optimal NWST value (for any set of 
terminals) is also no more than the optimal MWCDS value. Hence, if there is an a-approximation 
for MWDS(or equivalently WUDC) and a /3-approximation for NWST, then there is an q; + /3 factor 
approximation algorithm for MWCDS. 

Zou et al. [41] show that there exists a 2 . 5 / 9 -approximation for NWST if there exists a p- 
approximation for the classical edge-weighted minimum steiner tree problem. The current best 
ratio for minimum steiner tree is 1.39 [7]. Thus, there exists a 3.475-approximation for NWST. 
Combining with our PTAS for WUDC, we obtain the following improved result for MWCDS. 

Theorem 37. There exists a polynomial-time (4.475 -|- e)-approximation for MWCDS for any fixed 
constant e > 0 . 

9.2 Maximum Lifetime Coverage in UDG 

The maximum lifetime coverage problem (MLC) is a classical problem in wireless sensor networks: 
Given n targets ti,... ,tn and m sensors si, ...Sm, each covering a subset of targets, find a family of 
sensor cover Si,... ,Sp with time lengths ti, ...,in [0, 1 ], respectively, to maximize ti + ... + Tp 
subject to that the total active time of every sensor is at most 1. MLC is known to be NP-hard [ 8 ]. 
Berman et al. [4] reduced MLC to the minimum weight sensor cover (MSC) problem through Garg- 
Konemann technique [24]. In particular, they proved that if MSC has a p-approximation, then MLC 
has a (1 -L e)p-approximation for any e > 0. Ding et al. [14] noted that, if all sensors and targets 
lie in the Euclidean plane and all sensors have the same covering radius, any approximation result 
for WUDC can be converted to almost the same approximation for MLC. Hence, the current best 
known result for MLC is a 3.63-approximation. Using our PTAS, we obtained the first PTAS for 
MLC. 

Theorem 38. There exists a PTAS for MLC when all sensors and targets lie in the Euclidean 
plane and all sensors have the same covering radius. 
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Let us mention one more variant of MLC, called maximum lifetime connected coverage problem, 
studied by Du et al. [18]. The problem setting is the same as MLC, except that each sensor cover 
Si should induce a connected subgraph. They obtained a (7.105 + e)-approximation when the 
communication ratios Rc is no less two times the sensing radius Rs- Essentially, they showed that 
an a-approximation for WUDC and a /^-approximation for NWST imply an a -|- /3-approximation 
algorithm for the connected MLC problem. Using our PTAS, we can improve the approximation 
ratio to (4.475 -L e). 

10 Final Remarks 

Much of the technicality comes from the fact that the substructures interact each other in a com¬ 
plicated way and it is not easy to ensure a globally consistent order. The reader may wonder what 
if we choose more than two disks (but still a constant) in a small square, hoping that the uncovered 
regions become separated and more manageable. We have tried several other ways, like choosing a 
constant number of disks in the convex hull of the centers in a small square. However, these seemed 
to only complicate, not to simplify, the matter. 

We believe our result and insight are useful to tackle other problems involving unit disks or unit 
disk graphs. On the other hand, our approach strongly relies on the special properties of unit disks 
and does not seem to generalize to arbitrary disks with disparate radius. Obtaining a PTAS for 
the weighted disk cover problem with arbitrary disks is still a central open problem in this domain. 
An interesting intermediate step would be to consider the special case where the ratio between the 
longest radius and the shortest radius is bounded. 
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A Missing Proofs 

A.l Missing Proofs in Section 4 

Lemma 39. The central angle of any uncovered arc is less than vr. 

Proof. We only need prove the arc of a gadget is less than vr. As we union all gadgets and add more 
and more disks in TL, the central angle of an arc only becomes smaller. So, now we fix a square 
gadget Gg(r). Consider the substructure above line DgDt. See the right hand side of Figure 1 for 
an example. Suppose an arc a with endpoints M and M' on the boundary of Gg. The center of 
the arc is in the central area of Gg. If the central angle ^(a) > vr, its center should lie in the cap 
region bounded by a and the chord MM'. Without loss of generality, we suppose M is closer to 
the line DgDt than M'. Draw a auxiliary line at M which is parallel to DgDt. If the line does not 
intersect the central area, the center a locates below the line MM' (Hence, outside the cap region), 
thus the central angle is less than vr. If the auxiliary line intersects the boundary of central area at 
point N. Then, we can see that 

\ MN \ = V|MD,|2-x2 + \ DsDt \ - ^/\ NDt \^- x ^, 

where x is the vertical distance between point M and line DgDt. It is not difficult to see that 
|MA| > 1. It means that the center of the arc locates below the line MN (Otherwise, the distance 
from the center to M is larger than |AM|, which is larger than 1, rendering a contradiction). So, 
the central angle of any arc is less than vr. □ 

A.2 Missing Proofs in Section 7.1 

Lemma 24 Consider two non-empty small squares T, T'. Suppose the square gadgets in the 
two squares are Gg(r) = (Ds,Di) and Gg'(r') = (d;,d( ). The active region pairs (Ari, Ar2) and 
{Ar'i,Ar' 2 ) are associated with gadget Gg(r) and Gg'(r') respectively. (Ari, Ar 2 ) and (Ar'^,Ar 2 ) form 
a double-mixture. Then, the following statements hold: 

1. Their corresponding squares T, T' are adjacent; 

2. The core-central areas of Gg{T) and Gg'(r') overlap; 

3. The angle between DgDt and D'^D'^ is 0(e) 

4 . None of the two core-central areas can overlap with any small squares other than T and T'. 

Proof. Suppose the core-central areas of Gg and Gg' are do and respectively. Because (Ari, Ar 2 ) 
and (Ar^, Ar^) form a double-mixture. There exists at least one disk centered in T' which appears 
in both Ari and Ar 2 . Thus, the disk is centered in Co- It means that Co and C(, overlap. Hence, T 
and r' are adjacent. See Figure 13. 

It is easy to see that any core-central area can overlap with at most two squares. Since Co and 
C'o overlap, at least one of them overlaps with both F and F'. Without loss of generality, suppose 
Co overlap with both F and F'. We only need to prove C(, cannot overlap with other squares other 
than F and F'. 

Our proof needs a useful notion, called apex angle. Consider a square gadget. See the left one 
of Figure 1 (and we use the notations there). The line DsE and DgF are the tangent lines to the 
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Figure 13: The case that arcs in two different active regions are not order-separable. The two 
adjacent squares are T = and T' = A 2 A 3 A 4 A 5 . (D^, D^) is the square gadget in T and 

) is the square gadget in T'. The active regions Ari and Ar2 belong to the Gg(r), while the 
Ar'^ and Ar 2 belong to the Gg'(r'). 


boundary of core-central area at point Dg. We define apex angle to be ZEDgF. We notice that 
0^0 only depends on the distance \DsDt\. When the side length of a small square is e, is 0(e). 
In fact, even the core-central area is completely determined by Dg and Dt- Thus, we can generalize 
the concept core-central area by using any two overlapping disks in T-l. 

Now go back to our proof. See Figure 13. The segments DgDt and DgD[ cannot intersect 
because they belong to different small squares. Suppose the squares F and F' are A 1 A 2 A 5 AQ and 
A 2 A 3 A 4 A 3 , respectively. A 2 A 3 is the common side. W.l.o.g, suppose Dg is closer to line A 2 A 5 
than D[. If any disk centered in can appear in Ari but outside the disk D(, (i.e., outside the 
disks U Dt U D'^), the core-central areas defined by D'^Dt ® and should overlap nontrivial (not 
only touch at point Z)'). Thus, the angle ZDtD'^D^ should less than 0(e). It means angle between 
DgDt and D'^D't is at most 0(e). Moreover, angle between line D'^D't and A 3 A 4 is no less than 
I — 0(e) > 0(e). Thus, cannot intersect A 3 A 4 . Hence, Gl), cannot overlap other squares. □ 

Lemma 25 Suppose active region pairs (Ari(Mi), Ar 2 (M 2 )) and (Ar'i(Mi), Ar 2 (M 2 )) are associated 
with gadget Gg and Gg' respectively. If Ai and Mi are order-separable, then M 2 and A '2 are also 
order-separable. 

Proof. We prove the lemma by contradiction. We only consider the arcs which have siblings. 
Suppose Ml and M'l are order-separable but M 2 and M 2 are not order separable. W.l.o.g., assume 
a -< a',\/a £ Mi,Va' G M'l. Since M 2 and M 2 are not order-separable, there exist arcs 61,62 G M 2 
and 6 ' G M 2 such that 61 ^ 6 ' -< 62 . (If not, there exist arcs 6 G M 2 and 61^,62 G M 2 such that 
b'l <b < 63 .) Suppose 02 is the sibling of 62 , a' is the sibling of 6 '. Thus, a' -< 02 based on the same 
proof to Lemma 13. It yields a contradiction to the assumption that a -< a',ya G Mi, Va' G M'l. □ 

and Dt do not locate in the same square. This is core-central area in the generalized sense.. 
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Lemma 28 Consider the substructure St which contains an active region of square F. The sub¬ 
structure can be cut into at most three smaller substructures, by doing label-cut twice. At most 
one of them contains the active region. Moreover, this new substructure (if any) is bounded by the 
region Dom(r"^) associated with F. 

Proof. See Figure 8. We use the same notations defined in Section 7.1. The entire active region 
[Ui-DeCo ^ sub-region of (D(P, 2) — — D^) n . Thus, any disk centered 

in D((5s, 1) n D{Qt, 1) n H~^ can cover the active region. So, if there is a disk centered in D(Qs, 1) n 
we add it into Ti such that the entire active region is covered. If not, we prove 
we can cut St such that the new substructure contains the active region is bounded by the region 
(D(T), 1) — Ds — Dt) n H~^. The disks centered in ]R(F) n D(Qs, 1) — D((5s, 1) H D{Qt, 1)) cover point 
Qg. Thus, the arcs of these disks are order-separable to the arcs in active region. Similarly, the 
arcs of disks centered in M(F) n D{Qt, 1) — D(Qt, 1) n D((5s, 1) are also order-separable to the arcs in 
active region. Thus, we can add two disks in 77 to label-cut St into at most three new substructures. 
Obviously, there is only one of them containing arcs centered in F — D(Qs, 1) U D{Qt, 1). Hence, 
this new substructure is bounded in the region (D(i7,1) — — D*) n H^. □ 

A.3 Missing proofs in Section 7.2 

Lemma 30 In each of the above iterations, one substructure St(b,A) is cut into at most two new 
substructures. Any original arc in A cannot be cut into two pieces, each of which belongs to a 
different new substructure. 

Proof. Figure 9 gives an explanation about the change of substructures before and after the pro¬ 
cess. Suppose the envelope is {ai[Qs,Qi],... ,ak[Qk-i,QtW- After we add the disk D(ai) in 
H, the baselines of the two new substructures are bi = b[Qs,Pj+i] U aj+i[Pj+i, Qi] and b2 = 
ai+i[Qi+i, Pf^i] U b[P/^^,(5i]. Since Oj+i is an arc lying on the envelope, there does not exist an 
arc b with endpoints {A, B) such that A -< Pj+i -< Pl_^_l -< B. So the endpoints of any arc in Sti 
cannot be in M(St 2 ). Thus, any arc cannot be separated into two substructures. □ 

A.4 Missing proofs in Section 7.3 

Lemma 40. Consider two disks D„ and D; in squares F„, F; respectively. Suppose ^ Gg(Fu), D/ ^ 
Gg(F0 and and D; intersect at points A and B. Suppose the side length of square is /r. If A AD ^B 
(or ZADiB) is more than 2^/^, Gg(F„) and Gg(F;) should overlap. Moreover, if intersection A 
(resp. B) is in U(77), the two arcs which intersect at point A (resp. B) are in the same substructure. 

Proof. Suppose and D/ intersect at point A and B. G Gg(F„), D; G Gg(F/). Thus, D„, 
locate in the same square F^ and Dg locate in the same square. The line DuDi intersect disk 
at point E and intersect disk at point F. \DuE\ and \DiF\ are no less than 1 — ^ according 
to the triangle inequality. Suppose AAD^B is 9. Then the length \DuDi\ equals 2cos(2). Since 
9 > 2^/2fl, we have \DuE\ -\- \DiF\ > \DuDi\, i.e. and overlap. Figure 14 illustrates the 
situation. If A (or B) is in U(77), it is obvious that the two arcs which intersect at A (or B) can 
cover the same point on the baseline (i.e., the intersection point of and D^) which means they 
are in the same substructure. □ 

Based on the Lemma 40, we prove Lemma 32 as follows. 
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Figure 14: and D; intersect at point A and B. Suppose the side length of square is If the 

central angle of arc a[A, B] is more than 2-^/2^. D*, should overlap, where G Gg(ru), G 

Gg(ro. 

Lemma 32 We are given a substructure St(b,^). a ^ A and two arcs 61,62 ^ A. 4 / 61,62 cover 
the same point on a, 61,62 should belong to the same substructure. 

Proof. We prove the lemma by contradiction. Suppose the arcs 0 , 61,62 are parts of D, Di,D 2 
respectively. Consider a point P on D that can be covered by all of o, 61 , 62 . D and Di intersect at 
Ai and Pi, meanwhile D and D 2 intersect at A 2 and B 2 . See Figure 15. Since D and Di belong 
to different substructures, we know AAiDBi is less than 2^/2p based on Lemma 40. Similarly, 
ZA 2 DB 2 is less than 2y^2|U. It is easy to see that ZD 1 DD 2 is less than 2y^2/i. When ZD 1 DD 2 is 
0{y/e) when /i = 0(e). IP 1 P 2 I is less than 1 as \DDi\ and \DD 2 \ are more than 1 and less than 
2. Suppose ^43 and P 3 are the two intersection points of Di and D 2 . Thus, we have ZA^DiB^ > 
7 r /3 > 2^/2^. Moreover, P is in By Lemma 40, we can see that it contradicts the fact that 

Di and D 2 belong to two different substructures. □ 



Figure 15: Consider a point P on D. 0 , 61,62 which belong to disks D,Di,D 2 respectively, can 
cover the point P. D and Di intersect at Ai and Pi, meanwhile D and D 2 intersect at A 2 and P 2 . 
Then, Di and D 2 belong to the same substructure. 
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